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In a first order theory of dissipative hydrodynamics, we have simulated hydrodynamic 
evolution of QGP fluid with dissipation due to shear viscosity only. Simulation confirms 
that compared to an ideal fluid, energy density or temperature of a viscous fluid evolve 
slowly. Transverse expansion is also more in viscous fluid. We also study the effect of 
viscosity on particle production. Particle production is enhanced, more at large py. The 
elliptic flow on the other hand decreases and shows a tendency to saturate at large px- 



1. Introduction 

A large volume of experimental data from Au+Au collisions at RHIC are suc- 
cessfully analysed in an ideal hydrodynamic model H However, experimental data 
do show deviation from ideal behavior. The ideal fluid description works well in al- 
most central Au+Au collisions near mid-rapidity at top RHIC energy, but gradually 
breaks down in more peripheral collisions, at forward rapidity, or at lower collision 
energies, indicating the onset of dissipative effects. To describe such deviations 
from ideal fluid dynamics quantitatively, requires the numerical implementation of 
dissipative hydrodynamics. 

Though the theories of dissipative hydrodynamics has been known for more 
than 30 years , significant progress t oward its numerical implementation has only 
been made very recently I 4 | 5 | 6 | 7 | 8 [9] ^ Variable Energy Cyclotron Centre, 
Kolkata, we have developed a numerical code (AZHYDRO-KOLKATA) to solve, the 
1st order dissipative hydrodynamics in 2+1 dimension (assuming boost-invariancc 
in the longitudinal direction) and currently extending the code to 2nd order dissi- 
pative hydrodynamics. Some of the results from 1st order viscous hydrodynamics 
will be presented here. 



2. 1st order dissipative fluid dynamics 

Relativistic dissipative hydrodynamics has been discussed in detail in refl ^^ l i n 
dissipative hydrodynamics, dissipative fluxes are assumed to be small. The entropy 
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current is expanded in terms of dissipative fluxes. In the 1st order theory the 
expansion contain terms linear in dissipative fluxes, explicit form of which could 
be obtained by satisfying the second law of thermodynamics d^S^ > 0. 1st order 
theories are acausal, signal can travel faster than light. This is corrected in 2nd 
order theory entropy expansion contains terms quadratic in dissipative fluxes. 
Naturally, 2nd order theories are more complex. 

We consider a QGP fluid in the central rapidity region with net zero baryon 
density (tib = 0). We also neglect all the dissipative effects (e.g. heat conduction 
and bulk viscosity) other than the shear viscosity. We work in the Landau energy 
frame. Energy-momentum tensor, including the shear pressure tensor is written 
as. 

T"" = -pA^+Tt^" (1) 

where is the hydrodynamic 4-velocity (u^u^ = 1) and A^" = g^ v — u^u v is the 
projector, orthogonal to . e = u^T^Uu is the energy density, p is the hydrostatic 
pressure. T^ u satisfies the conservation law, 

g^u = q ( 2 ) 

In the first order theories , the shear stress tensors are written as, 



- (A^A^+A^A^) - -A^A^ 
2 y '3 



(3) 



where 77 is the shear viscosity coefficient. 7r MI/ is symmetric (ir^ v — 7T U,± ), traceless 
(ir^ = 0) and transverse to hydrodynamic velocity, (m /i 7t a " / = 0). The 16-component 
7r A " y has only 5 independent components. 

Heavy ion collisions are best described in terms of proper time r = \Jt 2 — z 2 
and rapidity r] s — ^lnfirf - In ( T TX,y,i] s ) coordinates, with longitudinal boost- 
invariance, the hydrodynamic 4-velocity can be written as, 

= (u T ,u x ,u v ,u Vs ) = ( 7± , 7x^,7x^,0), (4) 
with j± = l/ ^/l — v 2 — v 2 . The energy-momentum conservation equations are, 

d T T TT + d X (T TT V X ) + dy{T TT V y ) =~{P+ (5) 

d T f TX + 8 x (f TX v x ) + d y (f TX v y ) = -d x (p + n xx - n TX v x ) - d y (w xy - n TX v y l6) 
d T f T y + d x {f^v x ) + d y (f T y Vy ) = -d x (n x y - ^v x ) - d y (p + Jr™ - ^ Ty v y )(7) 

where v x = T TX /T TT and v y = T Ty /T TT , and we have used the notation "tilde" to 
represent quantities multiplied by the factor r, p — rp and similarly T y = tT % K 
We note that unlike in ideal fluid, in viscous fluid dynamics conservation equations 
contain additional pressure gradients containing the dissipative fluxes. Both T TX 
and T Ty components of energy-momentum tensors now evolve under the influence 
of additional pressure gradients. 
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With boost-invariance, number of independent shear stress tensor reduces to 3. 
tt tt , tt xx and n yy can be chosen as independent components (there could be other 
choices also^, 



T yy 



2 V 
2 V 
2ry 



-d x (^v x )--D( 1 lvl) + -(l + 1 lvl) 



(8) 
(9) 
(10) 



where 13 = = 7±(<9 T + v^c^ + v y d y ) is the convective time derivative, and 

is the local expansion rate, 6 = 2i + 9 r 7^ + d x (u x 7j_) + d y (v y ~/±). 

The other shear stress-tensors required in solving Eqs l5l6l7l can be obtained from 
the constraints satisfied by 7r A " / , (i) tt^ = 0:tracelessness, (ii) u^tt^ — 0: transverse 
to u". 



tV = 


7T 


— IT — 


TT yV 


2v x n XT = 


7T 




2 vy 
— V y 7T 


2v y ir yT = 


^ rr 
7T 


- v> xx 


+ v 2 y K yy 


2v x v y iT xv = 


„^TT 
7T 


- v> xx 


- v 2 y ir yy 



(11) 

(12) 

(13) 
(14) 

Shear stress-tensor components contains time derivatives of velocities v x and 
Vy. Thus at time step r, one needs the still unknown time derivatives. In 1st order 
theories, this problem is circumvented by calculating the time derivatives from the 
ideal equation of motion , Du^ = and De = — [e +p)V M u M . 



3. Equation of state, viscosity coefficient and initial conditions 

Through the equation of state, the macroscopic hydrodynamic models make contact 
with the microscopic world. We have used the equation of state, EOS-Q, developed 
in refffl It is a two-phase equation of state. The hadronic phase of EOS-Q is modeled 
as a non-interacting gas of hadronic resonance. The QGP phase is modeled as that 
of a non-interacting quarks (u,d and s) and gluons, confined by a bag pressure B. 
Adjusting the Bag pressure, the two phases are matched by Maxwell construction 
at the critical temperature, T c = 164MeV. 

Shear viscosity coefficient (rj) of dense nuclear (QGP or resonance hadron gas) 
is q uite uncertain. In perturbative regime, shear viscosity of a QGP is estimated 
M r] = 86.473^ j^=t- With entropy of QGP, s = 37f^T 3 and a s wQ.5, the ratio 
of viscosity over the entropy, in the perturbative regime is estimated as, 

(-) w 0.135, (15) 

V S / pert 
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However, QGP produced in nuclear collisions is non-pert urb atiy e. It is strongly 
interacting QGP. Recently, using the ADS/CFT correspondence El shear viscosity 
of a strongly coupled gauze theory, N=4 SUSY YM, has been evaluated, rj = ^N^T 3 
and the entropy is given by s = ^-N^T 3 . Thus in the strongly coupled field theory, 

= -L « 0.08, (16) 

S / ADS/CFT 47T 

To demonstrate the effect of viscosity on flow and subsequent particle produc- 
tion, we use both the perturbative and ADS/CFT estimate of viscosity. 

Solving Eqs l5l6l71 require initial conditions. In the present demonstrative cal- 
culations, we have used the similar initial conditions as in refffl We just mention 
that in re fffl initial transverse energy is parameterised geometrically, with 25% 
hard scattering. At initial time Tj=0.6 fm, the central entropy density is Si„i=110 
/to -3 . Fluid velocities are assumed to be zero initially. In dissipative hydrodynam- 
ics, dissipative fluxes need to be specified also. We assume that by the equilibration 
time r,-, the dissipative fluxes attained their longitudinal boost-invariant values. The 
independent components at initial time are, ir xx = 2rj/Ti, ir vy = 1r\jTi and 7r rr = 0. 




4. Evolution of the viscous fluid 

To demonstrate the effect of viscosity, with the same initial conditions, we have 
solved the energy-momentum conservation equations for ideal fluid and viscous 
fluid. 

In Fig[T] we have shown the constant energy density contour plot in x-y plane, 
after an evolution of 5.2 fm. The black lines are for ideal fluid evolution. The red 
and blue lines are for viscous fluid with ADS/CFT (r]/s— 0.08) and perturbative 
(77/,s=0.135) estimate of viscosity. With viscosity fluid cools slowly. Cooling gets 
slower as viscosity increases. To obtain an idea of transverse expansion of viscous 
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fluid, as opposed to ideal fluid, in FigJSJ we have shown the constant temperature 
contours in r — x plane , at a fixed value of y=0 fm. Transverse expansion is 
substantially enhanced in a viscous fluid. More the viscosity, more is the transverse 
expansion. The plot also indicate that at late time, fluid at x=y=0 behaves similarly 
to a ideal fluid. 

1st order theories are acausal, signal can travel faster than light. Acausality can 
lead to unphysical effects like reheating early in the collisions. However, we have 
not found any evidence of reheating. For small viscosity, with initial conditions as 
required in RHIC energy Au+Au collisions, effect of causality violation is minimum. 




012345 012345 
PT (GeV) p T (GeV) 



Fig. 3. Pt distribution of n for ideal and Fig. 4. Elliptic flow as a function of the trans- 
viscous fluid. verse momentum. 



5. Particle spectra and elliptic flow 

Viscosity generates entropy and particle production is enhanced. Viscosity influ- 
ences the particle production by (i) changing the freeze-out surface (freeze-out 
surface is extended) and (ii) by introducing a correction to the equilibrium dis- 
tribution function. For small departure from the equilibrium, the non-equilibrium 
distribution function can be approximated as, 

f(x,p)=fW(x,p)[l + <j>(x,p)], (17) 

(j)(x,p) is the deviation from equilibrium distribution function f 1 -- ^. With shear 
viscosity as the only dissipative force, (j)(x,p) can be locally approximated as, 

cj>{x,p) = Cn^p"; C = 2T2{ l +p) (18) 

completely specifying the non-equilibrium distribution function. The standard 
Cooper-Frye prescription for particle production from freeze-out surface (S^) can 
be employed to obtain the particle spectra, Detailed expressions for EdN/cPp with 
non-local equilibrium distribution function Eqs 1 171181 can be found in ref. 
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In FigO pr spectra for ir~ , from freeze-out surface at Tp = 158MeV is shown. 
Particle production is increased in viscous dynamics. We also note that effect of 
viscosity is more prominent at large pt than at low px- pr spectra of pions are 
flattened with viscosity. 

We have also calculated the elliptic flow in the model. Being a ratio, elliptic flow 
is very sensitive to the model. Experimentally, elliptic flow saturates at large pr- In 
FiglUpT dependence of elliptic flow is shown. Elliptic flow decreases with viscosity. 
As viscosity increases, elliptic flow reduces. We also note that both for ADS/CFT 
and perturbative estimate of viscosity, elliptic flow indicate saturation at large pr- 
The result is very encouraging, as experimentally also elliptic flow tends to saturate 
at large px ■ 

6. Summary and conclusions 

In a 1st order theory of dissipative hydrodynamics, we have studied the boost- 
invariant hydrodynamic evolution of QGP fluid with dissipation due to shear viscos- 
ity. In this model study, we have considered two values of viscosity, the ADS/CFT 
motivated value, rj/s «0.08 and perturbatively estimated viscosity, rj/s ssO.135. 
Both the ideal and viscous fluids are initialised similarly. Explicit simulation of 
ideal and viscous fluids confirms that energy density/temperature of a viscous fluid 
evolves slowly than its ideal counterpart. Transverse expansion is also more in vis- 
cous dynamics. For a similar freeze-out condition freeze-out surface is extended in 
viscous fluid. 

We have also studied the effect of viscosity on particle production. Due to vis- 
cosity particle production is enhanced, more at large px- The elliptic flow on the 
otherhand reduces and shows a tendency of saturation at large px- 
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